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Stellar populations of bulges at low redshift 


Patricia Sanchez-Blazquez 


Abstract This chapter summarizes our current understanding of the stellar popula¬ 
tion properties of bulges and outlines important future research directions. 


1 Introduction 


The stellar populations of bulges provide a fossil record of their formation and evo¬ 
lutionary history, including insights into the duration and efficiency of the primary 
epochs of star formation. In previous chapters, we have learnt that there are three 
main types of bulges: the classical, the disky, and the boxy/peanut. All of them show 
different structural and kinematical properties and different formation scenarios are 
proposed to explain them. In these scenarios, the star formation history is predicted 
to be different, e.g., the proposed mechanisms to form classical bulges imply rapid 
and efficient star formation while disky bulges are believed to form slowly and at 
low er redshifts from the inflow of mainly gaseous mate rial to the center of the galaxy 
(see lWvse et al.L Q 9971: iKormendv & Kennicuttl l2004 for reviews). 

Boxy/peanut bulges are believed to be parts of bars seen edge-on and have their 
origin in vertical instabilities of the disk. Therefore they are expected to have a sim¬ 
ilar stellar population compared to the inner disk. In principle then we should be 
able to distinguish between different formation scenarios simply by studying the 
different ages, metallicities and abundances ratios of the bulges. The situation, how¬ 
ever, it is not that straightforward; as discussed in other chapters (see e.g. Chapter 
2.3) internal processes related with disk instabilities can also occur at high redshift 
and in short timescales. Furthermore, bulges with properties resembling pseudob¬ 
ulges can form, not only by internal processes related with the pre sence of non- 
axisymmetric components, but also by accretion of gas and galaxies (IGuedes et al 
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2013t lOuereieta et al.L 12015t lEliche-Moral et all 1201 It lObreia et al.L 120131) . How 


ever, getting information on the ages, metallicities and abundances ratios in the dif¬ 
ferent types of bulges constitute, undoubtedly, a strong constrain for scenarios of 
bulge formation and, therefore, several authors have studied the problem. With the 
exception of the MW and M31, in which we can resolve individual stars, studies of 
bulges have to deal with integrated properties, through their mean color or absorp¬ 
tion lines. Such unresolved stellar populations studies have been far less common 
for bulges than for elliptical galaxies. The reason is that disk galaxies have more dust 
and ionized gas. The first affects the colors and the second fills the Balmer lines, the 
most important age diagnostics in the optical. In addition, bulges have, in general, 
lower surface brightness than ellipticals and the presence of several morphological 
components, such as disks, bars, rings, etc., complicates the interpretation of the 
results. Lastly, the light coming from the disk may contaminate the bulge spectrum 
in a way that is difficult to quantify. This problem is especially acute for studies of 
stellar population gradients. 

Furthermore, over many years, unresolved stellar populations studies have been 
done comparing the integrated colors or absorption lines with the theoretical predic¬ 
tions for single stellar populations (SSP); that is, an essentially coeval population of 
stars formed with a given initial mass function with the same chemical abundance 
pattern. While this scenario may not be a bad approximation for massive elliptical 
galaxies, bulges, especially those formed secularly, are believed to have a more ex¬ 
tended star formation history. This means that the young populations, which have 
low mass-to - light r atios, bias the analyses of composite populations, if present (e.g. 


Trager et al. . 200(i) . 


The relatively low number of studies, the small - and biased - samples, and the 
difficulties pointed out above have led to a lack of consensus about important results 
concerning the stellar populations of bulges, as I will show in this review. How¬ 
ever, in the last decade, stellar population models which predict not only individual 
spectral features , but the entire synthetic spectra for a population of a given age 
and metallicit^ Vazdekisl 199^ Bruzual & Chariot! 20031 Vazdekis et al. , 20 id 
Coelho et al.L 2007t Wale her et al. . 2009; Conroy et akf 2014 ) have been released. 
The availability of these models is stimulating the development of numerical algo¬ 
rithms to invert the observed galaxy spectrum onto a basis of independent compo¬ 
nents (combination of single stellar populations, age-metallicity relation, and dust 
extinction). Also, new specialized software allow the se paration of t he lig ht com¬ 
ing from the stars and ionized gas in a reliable way (e.g. Sarzi et all 20061). In ad- 


dition to this, new data from integral field spectrographs (e.g. iBacon et al.L 12001 


Cap pellari_&_e^l.l 201 ll Blanc et all 2013) a re cha nging the way we see galaxies 
dSanchez & et al.Ll2012HRosales-Ortega et al. . 2010l) . The analysis of these datasets 
allows one to associate stellar population properties with morphological and kine- 
matical characteristics of the galaxies, making the interpretation of stellar popula¬ 
tions more secure. Therefore, the development of the field is very promising and we 
foresee important advances in the decades to come. 

In this chapter, I will try to review the state of the art in the area, trying to high¬ 
light the necessary steps to get a better understanding of the star formation histories 
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of these complex systems. Section |2] summarizes the general results obtained with 
single apertures. Section [3 compile the works on the possible influence of bars in 
the stellar populations of bulges and on the stellar populations of bars themselves. 
Section|4]outline the results obtained with full spectral fitting techniques and Sect.|5] 
reports on the studies of stellar population gradients in bulges. In Sect. |6] I show 
the results about the possible connection between the stellar populations of bulges 
and disks while in Sect.|7]the main results are summarized. In Sect.[8l I give some 
thoughts of what I think the next steps for the study of stellar populations in bulges 
should be. 


2 Results obtained with single aperture 


2.1 General properties 


The first studies of stellar populations in bulges were performed using optical 
and near-infrared broadband photometry (Balcells & Peletierl 1994 ; Terndrur) et al 


1994 ; Peletier & Balcells , 19961 Bell & de Jong , 200C ; de Jong , 1996 ). These works 
demonstrated that changes in the bulge colors are linked to galaxy luminosity, po¬ 
tential well, and local surface brightness, with more massive/luminous bulges and 
higher surface brightness regions being redder than less massive/luminous and lower 
surface brightness ones. They also showed that early-type bulg es are red, as red as 
ellipti cal galaxies, and with very little dispersion in their colors (IPeletier & Balcellsl 
1996h . These results do not apply to the few late-type galaxies analyzed, where sig¬ 
nificantly bluer colors are measured. 

Early interpretation of these data pointed to early-type bulges being as old as 
ellipticals with late-type, less massive ones being younger and/or more metal poor. 
The small dispersion in the colors was interpreted as being due to small dispersion in 
the age of early-type bulges (maximum of '^2 Gyr). This, in principle, is in agree¬ 
ment with the classical and disky bulge formation scenarios (see Sect. [T]i if, as it 
seems the case, secula rly formed disky bulges are more common in late- type and in 
less massiye galaxies ( Kormendy & Kennicutd l^04HGanda et alll2009lpl 

One problem of using colors is the well known age-metallicity degeneracy 
( Worthey , 19941) . Bluer bulges can be either younger or more metal poor and, with¬ 
out this information, it is difficult to extract conclusions about their formation mech¬ 
anisms. Further complications are the presence of emission lines and dust extinction 
that also affect the colors. In particular, dust extinction depends on the inclination 
and, therefore, inclination is another parameter that needs to be taken into account 
when comparing the colors of different types of bulges and also when comparing 


* Note, tiowever, that a significant number of local massive spiral galaxies appears to have dominant 

i iseudobulges, defined a s those bulges with n < 2, that includes both disky and boxy/peanut bulges 
Kormendy et allllOlOl) . Furthermore, pseudobulges are also found in SO galaxies (Laurikainen et 
al. 2012). 
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the results from different studies (see Ganda et al. . 2009f) . All these obstacles make 
very difficult to extract useful conclusions about the stellar populations of bulges 
using only colors. For these reasons studies with colors need to be complemented 
with those using information of the absorption lines with diff erent sensitivities to 
age and metallicity and that are not affected by dust extinction ( MacArthur . 2005h . 
The first spectroscopic studies of bulg es analyzed the relation of line-strength in- 


dices (the so-called Lick/IDS indices, see lGorgas et al.[ri993HWorthev et aU 1 19941) 


with the central velocity dispersion (a hereafter) - used as a proxy for the dynam¬ 
ical mass of the galaxy. Lick/IDS indices measure the strength of the most promi¬ 
nent absorption lines in the optical galaxy spectra and are sensitive to changes 
of the mean age, chemical abundances and, to a lesser ex t ent, the initia l mass 
function (e.g. \4zdekis_etakL 201^ Bruz ual & Chariot , 2003 ; Thomas et al. , 2003 


SchiavonL 120071: IConrov & van Dokkumi l2012h . These studies confirmed the sim- 
ilarity of bulges w i th elliptical galaxies for early-type g alaxies (earlier than She, 


Bender et al.. 1993; Fisher et al.. 19961: Idiart et al.. 199^ 


However, it has been pointed out that this similarity may be due to the fact 
that the majority of these early analyses were pe rformed on samples that are bi - 
ased towards early-type spirals (earlier than She, Kormendv & Kennicutd 2004). 
In the last few years, however, several studies have included in their samples 
late-type bulges and analyzed, mostly, the relation between the Mg-sensitive in¬ 
dices (Mg 2 and Mgb) and the central <J. When these bulges were included, dif¬ 
ferences between elliptical and bulges were found. However, the nature of these 
differences is still not clear. Some authors claim that bulges are located below the 
Mg-a relation obtained for elliptical s, which is commonly interpreted as bulges 
having a younger s tellar populat i on (Prugniel et al.L 200 it Chiappini et al. . 2002; 


Ganda et al.L l2007t iMorelli e t al.L 120 081). Other authors find t hat the slope of the 


Mg -g is steeper for bu l ges (filcon-Barroso et ak . 2002), while Trager et al. ( 1999h 
and IProctor & SansomI (l2002h report that only low-mass bulges depart from the rela¬ 
tion between spectral indices and g drawn by large bulges. On the other hand, other 
studies do not find any systematic difference in the Mg-g relation of bulges and el¬ 
liptical galaxies but find that the scatter among this rela t ion for bulges is larger tha n 


liptical galaxies but tind mat tne scatter among mis rela t ion tor bulges is larg i 
the equivalent one for ellipticals ( Moorthv & Holtzman . 2006t Peletier et al.L 


20071). 


Similar conclusions were obtained using other line-strength indices. 

Some of the discrepancies in the conclusions of different studies may be due to 
differences in the mass distribution (or central g) of the selected sample. For exam¬ 
ple, Fig. [T] shows the relation between the Mgb and Hj3 indices (measured in mag- 
nitudes) for a samp le of late-type and early-type bulges and elliptical galaxies from 
Ganda et"^ (l2007l) . It can be seen that despite the relations between line-strength 
and indices followed by late-type bulges and early-type galaxies running apparently 
parallel to each other, the differences may just be related with the different range of 
central g. In fact, the differences disappear at low-g, where even early-type galaxies 
(E and SC0) deviate from the relation defined by massive ellipticals. 


^ Clas.sically, studies of stellar populations include E and SO galaxies in the same group. 
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log(a) (km/s) 
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Fig. 1 Line-strength indices Mgb' and Hji' expressed in magnitudes, agains t central velo c ity dis - 
persion. The black symbols represent the sample of late-type bulges from iGanda et alj i2007h . 
and the blue and the red symbols, the sample of Sa galaxies and E and SO respectively from the 
SAURON survey. T he dotted and dashed-do t ted bla ck lines over-plotted in both panels are the 
relations obtained bv ISanchez-Blazguez et alj i2006jl for low- and high-density environments, re¬ 
spectively, and the green solid lines are the relations detennined using the late-type galaxy sample. 
Representative error bars are added at the bottom right of each panel; the black one refers to both 
the e arly- and late-type spiral samples, while the red one refers to the E/SO galaxies. Figure taken 
from iGanda et alj iioOTh . 


Furthermore, bulges, contrary to massive elliptical galaxies, are rotationally sup¬ 


porte d. Some authors have cautioned (IPrugniel & SimienLll994tlFalc6n-Barroso et al 


20021) that by not taking into account the rotation in the cj measurements, one may 


be underestimating their binding energy. T he contribution to the rota tion may by 
calculated as 0.51og(l +0.62y^/g^), as in Pmgniel & Simien ( 19941) . where V is 
the rotational velocity. iFalcon-Barroso et"^ ( 2002h claim that a mean V jo ~ 0.5 
suffices to bring the bulges back to the Mg2-<7 relation defined by giant elliptical^. 
Other studies have also claimed a better correlation between the line-strength in¬ 
dices and Vniax (an indica tor of the total potent ial well) than between line-strength 


indices and the central a (IPrugniel et al. L 120011) . 


2.2 Comparison with SSP models 

Nevertheless, the similarity found by some authors between the index-cr relation of 
bulges and elliptical galaxies may not reflect a real similarity in their stellar content. 
Line-strength indices are not free from the age-metallicity degeneracy (e.g., Mg 2 
and Mgb can be lower in younger or in more metal stellar populations). Therefore, 
large differences between the ages of bulges and ellipticals could exist and not be 


^ Note, however, that the majori ty of low-luminosity elliptical galaxies are also rotationally sup¬ 
ported lEmsellem & et alil201 ih . 
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reflected in these re lations if there is a complementary age-metallicity relation (e.g. 
Trager et al. . 2000l) . The advantage of using these characteristics, though, is that 
the sensitivity to variations of age and metallicity of each different index varies. A 
way to partially break the age-metallicity degeneracy is to combine indices more 
sensitive to mean age variations (i.e., the Balmer lines) with those more sensitive to 
abundance variations in the so-called index-index diagrams. Figured shows one of 
these diagrams combining the composite index 

[MgFe]' = \/Mgb(0.28Fe5270 -1- 0.72Fe5335), 


which is fairly insensitive to variations of [a /Fe] abundance^ dThomas et al. , 2003 ) 
and the Balmer index Hj3. 

Several authors have used this technique to compare the index values with the 
predictions of SSP models. These comparis ons show that bulges have a large range 
in SS P-equivalent ages from ^2 to 13.5 Gyr ( Peletier et al. , 2007 ; Moorthv & Holtzman , 


2006h and metallicities. They also report a correlation between both the SSP- 


equivalent age and metallicity and central <7. In general, they found that more lumi¬ 
nous/massive bulges were older and more metal rich. They also inferred that more 
massive bulges have a larger ratio of a-element^ with respect to Fe, which is usu- 


The relations were similar to those found for elliptical galaxies ( 

Bica. 19881 Jablonkaet al.. 

19961 

Idiart et al.l 19961 Casuso et al.l 19961 Goudfrooii et al. 

J1999I: Trager etal.L 

I999I 

Thomas & Davies. 20061 Moorthv & Holtzmanl 20061; Jablonkaet al., 20071: 

Ganda et al.. 20071. On the other hand. Proctor & Sansoml ( 2002h and Prugniel et alJ 


(l200lh found that, contrary to what happens in elliptical galaxies, both Fe and Mg 
were correlated with a in bulges, resulting in the lack of a tight correlation between 
Mg/Fe (a proxy for [a/Fe]) and a. This result needs to be corroborated by other 
studies. 

Similar to the results obtained with line-strength indices, the comparison of SSP- 
equivalent parameters of bulges and ellipticals reveal that both have very similar 
properties, at least in samples of bulges earlier than She. Figure [3 illustrates the re¬ 
lation between the SSP-equivalent ages, metallicities, and [g/Fe], and the central 
(7 for bulges, SO, and elliptical galaxies ( Thomas & Daviea l2006h . It can be seen 


that, at a given central (7, the stellar population parameters of bulges and elliptical 
galaxies are indistinguishable. These results support the idea that bulges (with mor- 


^ The a elements are those chemical elements predominantly formed via fusion with a helium 
nucleus. Their most abundant isotopes therefore have nucleon numbers that are multiples of four 
(e.g., O, Ne, Mg, Si, S, Ar, Ca, Ti). These elements are mainly synthesized in Type II supemovae, 
while Type la supemovae produce elements of the iron peak (V, Cr, Mn, Fe, Co and Ni). In chem¬ 
ical evolution models, type II supemovae produce an early enrichment of a-elements followed by 
a subsequent enrichment of iron-peak, Type la supemovae products. In the absence of other mod¬ 
ifying factors, this implies that [cc/Fe] can be used as a ’galactic clock’ for the duration of the star 
formation. 

^ What it is usually measured is the Mg abunda nce through the Mg b index. Other g-elements , 
like Ca or Ti may f ollow different patterns (e.g. IConrov et all 120141 : [Graves & SchiavoTl l2008l : 
ICenarro et alll200^ . 
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[MgFe]' 


Fig. 2 H6 vs. [MgFe] ' in the central regions of bulges and ellipticals. Magent a crosses are galax¬ 
ies fro m IXrager et^y jjj99|f Green symbols ai’e bulges and ell ipticals from IProctor & Sanson] 
i2002h . The MW jPuzia et alil200 2h and M3 1 jPuzia et al.L[200^ are shown as ’+’ symbols. Blue 
and red symbols are from iMoorthv & HoltzmanH2006ll . The color of the symbol is chosen accord¬ 
ing to whether they are redder or bluer than B-K=4. Larger sizes indicate larger central velocity 
dispersion. Bu lges for which there is no color infonnation are in black. Solid lines represent the 
predictions of iThomas et al.l i2003h for SSPs of constant metallicity (as indicated in the labels) 
while dotted lines represent the predictions for populations of constant ag e, with age increasing 
toward s the bottom of the panel. For other details regarding the hgure see iMoorthv & HoltzmanI 

<2006ll . 


phological types earlier than She) were formed with very little influence from the 
disk, in a process similar to the one that formed elliptical galaxies. 


In s ampl es that contained late-type gal axies (late r than She), bothiPrugniel et al 
( 200 ih and Moorthv & HoltzmanI ( 20061 see also iMorelli et al.l l2008li And three 
types of bulges in the comparison of line-strength indices and stellar population 
models: the old-metal rich (OMR), the young metal rich (YMR) - which are bulges 
with ages less than 3 Gyr and super-solar metallicities - and the metal poor (MP), 
with sub-solar metallicity. These classes seem to be sensitive to the Hubble type. 
All the early-type (SO-Sab) bulges are metal-rich. The red early-type bulges are in 
the OMR region while the blue early-types reside in the YMR region. Metal-poor 
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Fig. 3 Stellar population parameters ve rsus central velocity dispersion. Open circles are early-type 
galaxies from iThomas & DaviesI i2006h . ellipses are spiral bulges with ellipticity increasing for the 
later-types (see labels in the right-hand bottom panel) and the filled s quare is the integrate d light of 
the MW bulge. Small grey-filled circles are early -type galaxies from ^omas et al] l l200.'5ll . Central 
stellar populations are shown. Figure taken from iThomas & DaviesI ■ 2006h . 


bulges are all late-types, but late-type bulges are found in all three regions. A com¬ 
parison of the SSP properties of late-type bulges and ellipticals galaxies at a similar 
(7, however, remains to be done and, therefore, as it was the case with the line- 
strength indices, it is not clear if late-type bulges are younger, for being late-type, or 
for having low a. 


2.3 Relation between stellar population and structural properties 


A more direct way to test the different proposed scenarios for bulge formation is to 
compare the stellar population properties of a sample of bulges with morphological 
or dynamical properties that distinguish them as classical, disky, and boxy/peanut. 
As we have seen in previous chapters (see Chapter 2.2), there are several observables 
used to separate bulges and pseudob ulges and different a uthors have employed dif¬ 
ferent properties to perform this task. Carollo et al. ( 200lh analyzed V, H, and J HST 
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images of a sample of bulges with exponential (typical of disky and boxy/peanut 
bulges) and R*/'* luminosity profiles (typical of classical bulges), finding the former, 
on average, bluer than the latter (by ~0.4 mag in <V —H >), which could be the 
consequence of a younger and/or of a lo wer metallicity stellar popu lation. They also 
found, in agreement with the results of Peletier & Balcell^ (Il996 ). that the colors 
of those bulges showing a profile were red and very homogeneous, while for 
the exponential bulges the scatter was significantly largeiQ. They interpreted these 
results as a delayed formation of the exponential bulges compared with those hav¬ 
ing an R*/^ profile, which formed their stars in the early Universe. IProrv & Fishei 
( 2007h used a different approach and separated classical and pseudobulges mor¬ 
phologically. Pseudobulges were those showing nuclear bars, nuclear spirals, and/or 
nuclear rings and classical bulges those featureless structures rounder than the outer 
disk. Separating the bulges in this way and comparing with their visual morphologi¬ 
cal types, they studied the location of bulges in the color-magnitude diagram, finding 
that Sc galaxies and later types do not contain classical bulges and are located al¬ 
most entirely in the blue cloud in the color-magnitude diagram. Intermediate Sa-Sbc 
type galaxies, on the contrary, contain both classical and pseudobulges. While 87% 
of the galaxies with pseudobulges were in the blue cloud, all galaxies with classical 
bulges were in the red-sequence. These authors stress that the differences in colors 
are not due to a different contribution of bulge and disk to the total galaxy color as 
blue and red galaxies share a range in bulge-to-total ratio. This is shown in Fig. 01 
where the global color of the galaxies is plotted as a function of the bulge-to-total 
(B/T) ratio. It can be seen that both bulges and pseudobulges coexist in the region 
of B/T values ranging from 0.05 to 0.45 and that, in this region, classical bulges are 
redder than pseudobulges. They also found that, in general, pseudobulges classified 
with the morphological features described above were more diffuse and had lower 
Sersic indices than classical bulges. 

Another wa y of separating classical bulges from pseudobulges was adopted by 
Gadottil ( 2009h who ident ified pseudobulges as those lying below the Kormendy 


relation ( Kormendvl 1977 ) defined by elliptical galaxies. He established that pseu¬ 
dobulges defined this way were, in gen eral, 0.2 mag bluer in the (g — i) color than 
the classical bulges. On the other hand, Fernandez Lorenzo et al. (12014 ). who used 
the Sersic index to differentiate between classical and pseudobulges, found the latter 
as red as the former at the same luminosity. Only in fainter pseudobulges they did 
measure bluer colors. 

Again, colors are affected by the age-metallicity degeneracy and dust extinction 
and, in principle, more information may be obtained from the study of absorption 
lines and their comparison with stellar population models. Several works have also 
com pared the line-s t rength indices of bul ges with different structural characteris¬ 
tics ( Williams et al. . 2012 : Gadottil 20091) . The first studied a sample of edge-on 
boxy/peanut bulges findi ng that they follow the same central index-a relation as el¬ 


liptical galaxies (see also lJablonka et all 120070 . altho ugh the s ample was biased to¬ 


wards early-type galaxies (SO-Sb). On the other hand, Gadotti ( 20091) found a strong 


® Note, however, that many of the exponential bulges were showing colors as red as the R*/^ bulges. 
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Fig. 4 Distribution of bulge-to-total ratios, B/T, of intermediate type (Sa-Sbc) galaxies with pseu¬ 
dobulges (blue triangles) and classical bulges (red circles) with respect to their glob al u-r color. The 
dashed line marks B/T=0.45 above which only classical bulges are found. From IProrv & Fisheil 
j2007ll . 


correlation between the D4000 index and the Sersic index (n), indicating younger 
populati on^ in galaxies with to wer n. On the contrary, for a sample of early-type 
galaxies, IVazdekis et al.l (120041) did not hnd a correlation between n and age, while 
they found a strong correlation between n and [Mg/Fe]. The different behavior for 
bulges and ellipticals is very interesting, but more studies of this kind using larger 
samples are still needed to conhrm or refute the t rends. Note tha t the on ly study tar¬ 
geting specihcally boxy/peanut bulges is that of IWilliams et al.l (l2012h . The rest of 
studies cited above usually include both disky and boxy/peanut bulges in the same 
category, the ‘pseudobulge’ sample. 

The current lack of consensus between studies may be due to different criteria 
to separate classical and pseudobulges. Some might include disky and boxy/peanut 
bulges in the same category, pseudobulge, without making any distinction between 
them. Furthermore, different distributions of the galaxy luminosities can also lead 
to discordant results. It is clear that low-luminosity, low-mass bulges are bluer than 
more massive and brighter ones, but it is not clear if, at the same luminosity, bulges 
with different structural characteristics share the same color. 


^ Actually the author did not compare the index with stellar population models and, therefore, a 
variation in metallicity is also possible. 
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Thanks to th e 2-dim ensional data of the SAURON survey (IBacon et all 1200Ih . 
Peletier et al. ( 2007h (see also Sil’chenko & Afanasiev . 2004ll noticed that, when 
present, young stellar populations in their sample of early-type bulges w ere concen¬ 


trated near the center , in dis ks or in annuli suggestive of resonance rings (IBvrd et al 


1 994) . iPele tier et al.l (l2007h realized that the studies comparing the line-strength in¬ 
dices of bulges and elliptical galaxies could be divided into two categories: those 
targeting inclined galaxies, which do not find any difference between the index-cj 
relation of bulges and ellipticals, and those sampling almost face-on galaxies, which 
find younger stellar populations in bulges compared with those of elliptical galaxies 
and a large scatter in the line-strength indices at a given a. The differences are espe¬ 
cially visible in galaxies with low a. They argue that bulges are composite systems, 
with two or more types of bulges coexisting in the same galaxies. The classical 
bulge is composed mainly of an old and metal rich population and the disky and 
boxy/peanut bulge can be younger and contain more metal-poor stars (although it 
can also be old). The discrepant results obtained in samples of different inclinations 
can be explained, according to these authors, by the different contribution to the 
bulge light of different subcomponents (classical, disky, and boxy/peanut). If the 
young component is a disk, then it is concentrated in a plane and it would not be ob¬ 
served in edge-on galaxies. These young components, however, do contribute to the 
integrated light of less inclined samples (if they are limited to the central regions). 
The result is supported by the observation of central dips in the velocity dispersion 
maps in 50% of the galaxies of their sample. 

The coexistence of 2 or more types of bulges (classical, di sky, and boxy/peanut) 
in some galaxies has been pointed out by several authors (lAthanassoulal 2005 


Gadottii 20091: Nowak et al. . 2010t Kormendv & Barentinei ioiOtlErwinT 2008h an d 

supported by theo retical studies dObreia et al.lbOlsHSamland & Gerhardl 2003 ). 


Obreia et'di ( 2013 ) propose a picture were the centers of most early-type spirals 


contain multiple kinematic components: an old and slowly rotating elliptical-like 
component, and one or more disk-like, rotationally supported components which 
are typically young but can also be old. 

This ‘two component model’ also explains the properties of our MW bulge. 
The MW is considered to have a boxy bulge, yet increased evidence of an old, 
a-enriched stellar pop ulation that formed on a sho rt time-scale has resulted in a two 
component model (e.g. Tsupmoto & Bekkil 2012h . It has been sh own that two stellar 


1994 

; Feltzing & Gilmore 

. 2000l van Loon et al.l 20031 Groenewegen & Blommaert 

2005 

; IZoccali et al., 2006 

; iFulbright et all 2007: IZoccali et all 20081) and that the 


I2007h . 


even if age determinations through color- magnitude diagram shows that most bulge 


stars in the Galaxy are older than 10 Gyr (Ortolan! et al.L Il995t iFeltzing & Gilmore , 


200(i Zoccali et al. . 2006t Clarkson et al. . l2008l) . 

The task of isolating the stellar population properties of the different subcompo¬ 
nents forming a bulge is difficult. Still, it has been tried by some authors. For exam- 
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pie, Williams et al, ( 201 ih study the stellar populations of two edge-on boxy/peanut 
shaped bulges. They place the slit along the major axis and observe with three offset 
in parallel positions. They found that NGC 1381 has a boxy bulge, with stellar ro¬ 
tation neither cylindrical (as would be expected for bars seen edge-on) nor strongly 
non-cylindrical and with a double hump on the rotation curve. The galaxy shows 
a metallicity gradient but no age gradient and a positive [a/Fe]. They explain the 
properties of these galaxies in an scenario where NGC 1381 has the three classes of 
bulges. The classical bulge formed their stars rapidly and explain the general trend 
in [a/Fe] as a function of height, as disk light (with its lower [a/Fe]) contributes 
less and less to the integrated spectrum. The boxy appearance is explained by the 
simultaneous presence of a bar (which appears boxy in projection), and the double 
hu mp of the rotation curve hints at the presence of a small disky pseudobulge (see 


also Sil’chenko et al. , 2010h . 


3 The influence of bars in building up the bulge 


It seems clear that some bulges have central disks (iPeletier et all I2007h . often (but 


not always) with young stars, which is usually linked to disk gas inflow and cen¬ 
tral star formation caused by internal secular processes related with the presence 
of a b ar ( Friedli & Benz . 1995 : Norman et al. . 19961 : Noguchi! 2000! Immeli et al 


2004li . However, this central rotationally younger component does not necessar¬ 
ily form due to internal processes. Major and minor mergers and external accre¬ 


tion of gas may resul t in the formation of a disky bulge (e.g. iGuedes et al.L 12013 


Ouereieta et al.Ll2015h. This idea may be supported by some observational studies. 


For example, iKannappan et al.l (120041) found, in a sample of disky bulges selected 
to be blue and, therefore, with central young stars, that all of them showed signs of 
recent interactions. 

Hydrodynamical cosmological simulations predict that both secular and exter¬ 
nal processes contribute to create disky pseudobulges with similar characteristics in 
both cas es, rotationa l ly supported and wi t h a young and metal poor stellar pop¬ 
ulation ( Obreja et al. . 2013! Guedes et al. . 20131) . Eliche-Moral et all ( 201 ih also 
analyze the effects of minor mergers on the inner part of disk galaxies, finding 
this process to be efficient in forming rotational ly supported stellar inner compo¬ 


nents, i.e., disks, rings or spiral patterns (see a lso iDommguez-Tenreiro et al. 
Aguerri et ah . 200 1! Scannapieco et al.! 2010f) . 


1998 


A way to quantify the importance of secularly formed disky bulges is to compare 
properties of galaxies with and without bars. While it has been found that emis¬ 
sion is enha nced in early-type spirals with bars with respect to those early-type non 
barred (e.g. Ho et al.! 199?! Huang et al. , 1996! Alonso-Herrero & Knapen! 


Jogee et al 


2001 


2005! lEllison et al.! 1201 Ih . the eyidence supporting the bulge building 


by bars from the ages of its stars has proyen to be elusiye. Seyeral authors haye tried 
the comparison using samples of face-on galaxies, where it is easy to morpholog¬ 
ically identify the bar. The differences, howeyer, haye not been firmly established. 
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Moorthv & Holtzman ( 2006 ) and Perez & Sanchez-Blazaued ( 201 ih found hints of 


lower ages and higher metallicities in barred galaxies compared to their counter¬ 
parts in unbarred at a given 00. They also found higher [a /Fe] abundances in barred 
galaxies with central velocity dispersion 2.2 > log a (km/s) > 2.35, but the oppo¬ 
site for 2 > log a (km/s) > 2.2. At fixed a and Vmax, barred galaxies appear to have 
larger central values of [MgFe]0 (which can be used as an indicator of metallicity 
independent of [a/Fe], see above) than non-barred galaxies (or galaxies with ellip¬ 
tical shape bulges) of the same o or Fmax. The differences, however, were not very 
significant in a statisti cal sense. 

On the other hand. ljablonka et alJ OOOTh found no difference between the stellar 
population properties of edge-on barred and non- barred galaxies. However, it may 
be difficult to detect a bar in an edg e-on galaxy. Ide Lorenzo-Caceres et al.l (120 12h 
and ide Lorenzo-Caceres et akl (120131) analyzed the stellar populations in the center 
of double-barred early-type SOs and spirals finding some signs of gaseous flows and 
young stellar populations. This population was not very prominent though. Never- 
th eless, all the above studies were affected by poor number statistics. 

Coelho & Gadottil ( 201 ih observed 575 face-on bulges in disk galaxies, of which 
251 contain bars. They found that, for bulges with masses between lO'^^'^M© and 
1O1O.85 m 0, the distribution of ages in barred galaxies is bimodal with peaks at 4.7 
and 10.4 Gyr. This bimodality is not seen in non-barred galaxies of a similar bulge 
mass range. The age distribution of barred and non-barred galaxies is, per contra, 
similar for bulges of masses lower than IO'^^'^Mq (i.e., the differences are only seen 
in massive bulges). These authors did not find any difference in the metallicity dis¬ 
tribution of barred and non-barred galaxies. These results are summarized in Fig.|5] 
where the distribution of ages for several mass intervals is shown for samples of 
barred and non-barred galaxies. 

Therefore, it is still not clear if there are differences between the stellar pop¬ 
ulations of galaxies with and without bars. What seems clear though is that, if 
there are differences, they are only visible in massive, early-type galaxies, a re- 
sult that is supported b y other studies analyzing the molecular gas concentration 
( Sakamoto et al.L 19991) . This is often attributed to the fact that early-type gala xies 


host larger and stronger b ars than late-type galaxies (e.g.jBufa_&Comte 


l9m al¬ 


though, as pointed out in iLaurikainen et akl (2007) and lLaurikainen etal.l (l2004l) . a 
longer bar does not implies a stronger bar and, in fact, in early-type galaxies the bar 
induced tangential forces are weaker because they are diluted by the more massive 
bulge0. 

This does not mean, necessarily, that bars are not efficient agents in building 
up bulges. It is still not clear if bars are long-lasting structures or not. If bars are 


* Although the differences found by iMoorthv & HoltzmanI l2006l) in the Hj 8 - 0 ' relation between 
barred and non-baired galaxies disappear when Vmax is used instead of O'. 

^ This index is defined as [MgFe]' = ^Mgb x (0.72 Fe5270 -t- 0.28 Fe5335 in iThomas et all 
120031) . 

Alt hough the majority of authors does not distinguish different types of SO. ILaurikainen et all 
120071) also show that early-type SOs have shorter bars than later type SOs, i.e., the trend of longer 
bars for early-type morphologies reverse in this morphological subclass. 
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log Mj, 


71 w/o bars 
63 with bars 


log Mb, 


77 w/o bars 
59 with bars 






Fig. 5 Normalized distributions for bulge ages, for several mass intervals, as indicated. Distribu¬ 
tions for b arred and non-barred gal axies are shown in red and black lines, respectively. (Figure 
taken from ICoelho & Gadottiil20llh . 


not long-lasting structures but recurrent patterns ( Bournaud & Combesl 2002 ) then 
the fact that we do not find differences between barred and non-barred galaxies 
would not necessarily imply that bars are not important for secular evolution but, 
simply, that non-barred galaxies could have been barred in the recent past. How¬ 
ever, most numerical simulations show that, once formed^ bars are robust struc¬ 
tures ( ShetL&Sellwo^ 2004; Athanassoula eLaL, 2005; pebatfeta eLaL , 2QQ6; 


Berentzen et al. . 2007 : Villa-Vargas et al. . 2010l: Kraliic et al. . 2012t Athanassoula et al 

2013h . Furthermore, at least in massive disk galaxies, bars have the same stellar 
populat ion properties of bulges (old, metal rich, and fg/ Fe] -enhanced stellar popu¬ 
lations; |sinchez^Bliz£uezeraD|20ni|pirezetd]^^ which , in many cases, are 
very different from that of the disk (see Sect.^. This result also supports (although 
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it does not prove, see Sect. [3T‘h the idea that bars formed long ago. The longevity 


of bar s is also suggested in studies of the bar fraction evolution (e.g. lSheth & et al 


20081) . which find a similar bar fraction at z~0.8 to that seen at the present-day for 
galaxies with stellar masses M* >10** M©. In addition, non -axisymmetric struc¬ 
tures. such as nuclear spirals, can drive gaseous inflows (e.g. Kormendv & Fishen 


200-51) . which could dilute the differences between barred and non-barred galaxies. 


3.1 Stellar population of bars 


As the debate of the durability of bars is still open and its influence may be crucial 
for the formation of bulges, it is also important to study the stellar populations hosted 
by bars. Very few works, how ever, have dealt with this problem. 


Gadotti & de Souza ( 2006 ) obtained the color and color gradients in the bar re¬ 


gion of a sample of 18 barred galaxies. They interpreted the color differences as 
differences in stellar ages and conclude t hat younger bars were hosted by galax¬ 
ies of later types (see also Gadotti . 2008h . However, as we mentioned in Sect. |2] 
the effects of age-metallicity degeneracy and dust extinction are strongly degen- 
erate in colors an d, th erefore, con c lusion s based on only colors remain uncertain. 
Perez et al. (2007) and Perez et al. ( 20091) performed an analysis of the stellar pop¬ 
ulation of bars in early-type galaxies using line-strength indices. They found that 
the mean bar values of SSP-equivalent age, metallicity, and [a/Fe], correlate with 
central a in a similar way to that of bulges, pointing to an intimate evolution of 
both components. Galaxies with high central <7 (>170 kms^*) host bars with old 
stars while galaxies with lower central velocity dispersion show stars with a large 
dispersion in their ages. 

These authors also analyzed the stellar population gradients along the bars and 
found three different behaviors; (1) bars with negative metallicity gradients. These 
bars have young/intermediate stellar populations (SSP-equivalent values < 2 Gyr) 
and have amongst the lowest stellar velocity dispersions of the sample; (2) bars with 
no metallicity gradients. These galaxies have, however, positive age gradients and 
(3) bars with a mean old stellar population and positive metallicity gradients (more 
metal-rich at the bar ends). 

The fact that bars are composed of old stellar populations does not mean that 
they formed long ago, as the bar might have formed recently out of old stars in the 
disk. One way to disentangle these two options is to compare the stellar popula tions 
of the disk and the bar at the same distances. In ISanchez-Blazauez et al.l(l201 ih . this 
comparison is made for two galaxies, finding that stars in the bar are older and more 
metal rich than those of the disk. Furthermore, the gradient in both parameters is 
much flatter in the bar. In general, they found that the stellar content of the bar is 
more similar to that of the bulge than to the disk. However, the sample of this study 


* * The stars in the bar can form in the disk long time ago, even if the bar have been recently formed. 
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remains small and biased towards early-type bulges. Clearly, a larger study sampling 
larger samples of galaxies covering all morphological types is still needed. 


4 Star formation histories 


The majority of stellar population studies in bulges are based on the comparison 
of the observed colors or spectral properties with the predictions of SSP models. 
While the SSP assumption may not be a bad approximation for massive elliptical 
galaxies, it is most likely not a good one for spirals, which are believed to have a 


more extended star formation history ( Kennicutt , 1993 : James et al. . 2008h . 


In cases were the star formation history has been more complicated than just 
a single burst, the interpretation of the results based on analyses of single-stellar 
population is difficult. Bulges with an intermediate SSP-equivalent age could have 
formed all their stars at intermediate epochs or almost all at very early times and a 
small fraction at recent epochs. 

To avoid these difficulties, some authors have tried to analyse their data as¬ 
suming more realistic star formation histories. There are several examples in the 
literature where t he stellar population properties have been derived using a para¬ 
metri c approach ( Ganda et al. . 2007 : MacArthur et al.L 2004; Kauffmann & et al 


2003h where a predefined shape for the star formation history and chemical enrich¬ 
ment is assumed. In this case, some parameters are fixed while others are fi t ted by 
comparing the observations with the predictions of the models. iGanda et al.l OOOTh 
compared three line-strength indices with the predictions based on two burst, and 
on exponentially declining star formation histories. They obtained, respectively, the 
age, metallicity and mass fraction of the youngest burst, and the e-folding (t) time in 
the case of the exponentially declining star formation history. In the first case, they 
concluded that degeneracies in the parameter space prevented them from extracting 
useful conclusions. In the second case, they found that bulges with a larger cen¬ 
tral (7 showed shorter e-folding times, more consistent with an instantaneous burst 
scenario. Low-a galaxies have larger T, indicating a more extended star formation 
history (see Fig. IS). 

The problem with this approach is that the results depend strongly on the priors 
- i.e., the wrong answer can be obtained with the wrong assumptions about the star 
formation history. Furthermore, the use of only a few line-strength indices makes it 
difficult to break the existing degeneracies in the parameter space, such as the age 
of the burst versus its strength, or the T-metallicity degeneracy. 

The avail ability of high-quality stellar libraries and associated stellar population 
mode l s (e.g. Sanchez-Blazquez et al., l2006bHVazdekis et al.Ll2O10HBruzual & Charloi 


2003 : Conroy & yan DokkumI 201^ that predict, not only indiyidual absorption 
line features but the whole spectral energy distribution, has allowed the deyelopment 
of new techniques that, fitting the whole spectrum, are able to obtain, not only SSP- 
equiyalent ages and metallicities, but a more, in principle realistic, star formation 
history. Furthermore, by considering the information proyided by the entire spec- 
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Fig. 6 Central aperture values for the e-folding time-scale t against central velocity dispersion a 
(both in units of decimal logarithm), in an exponentially declining star for mation scenar i o; the T 
values are obtained selecting amongst models with age=10 Gyr. Figure from iGanda et alJi2007h . 


trum, the age and the metallicity are more easily separated (ISanchez-Blazauez et al 


201 11) . These techniques are non-parametric - i.e., no predefined shape for the 


star formation history is assumed. Codes that have been used to s tudy t he stellar 

population propertie s of galaxies i nclude MOPE D ( Heavens et al.L |200(^: VESPA _ 

iToieiro et alll2007h: STECKMAP dOcvirk et alll2006bh: STARLIGH T dCid Fernandes et al 


2005h : SEDFIT dWalcher etZll2006ll and ULvSS dKoleva et al.ll200^ . Using these 
new tools, one can fit an observed spectrum in terms of a model built by a linear 
combination of a number of SSPs with different ages and metallicities. The kine¬ 
matics can be calculated at the same time by convolving the model with a Gaussian 
line-of-sight velocity distribution. In some cases, dust can be modeled assuming a 
reddening law. Figure Q shows an exa mple of a fit to the inte grated spectrum of 
NGC 628 using the code STECKMAP ( Ocvirk et al. . 2006bla ). together with the 
derived flux and mass fractions for sta rs of different ages and the age-metallicity 


relation (ISanchez-Blazauez et al.L 120141) . 


The problem of inverting a spectrum to derive detailed star formation and chem¬ 
ical enrichment histories is ill-conditioned - i.e., small fluctuations in the data can 
produce strong variations in the final solution - and different codes try to overcome 
these issue with different methods. The accuracy of the recovered star formation 
history depends critically on the signal-to-noise of the input spectra and, depending 
on this value, one can recover more or less different stellar populations described 
by an age and a metallicity. However, with spectra of enough quality, the meth¬ 
ods have shown that they can recover reliably both the age distribution and the 
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Fig. 7 Top panel: Example of the fit (red line) obtained with STECKMAP for the integrated spec¬ 
trum of NGC 628 (black line). The residuals from the fit (observed-fitted) are also showed. Hashed 
regions indicate those zones that were masked during the fit. Bottom panels: Mass, flux fractions 
and t he age-metallicity relation deriv ed with STECKMAP for the integrated spectrum of NGC 628. 
From lSanchez-Blazguez et^ l l2014h . 


age-metallicity relation (see, e .g. Cid Fernandes et al. . 2005 : Ocvirk et al. . 2006bl : 


Sanchez-Blazquez et alil20111) . 


Robust quantities, even when derived fr om a spectrum with l ow sig nal-to-noise. 


are the mean values of age and metallicity (ICid Fernandes et all 120051) that are usu 


ally weighted with the light (LW) or the mass (MW) of the stars. These are obtained 


< A >LW= 


If A; X flux,- 
I, flux, 


< A >MW = 


if A,' X mass,' 
I,'mass,' 


( 1 ) 

( 2 ) 
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log (Gyf) log <t>i,ght (Clyr) 


Fig. 8 Comparison of the SSP-equivalent age and metallicity as a function of mass-weighted 
(left-hand panels) and V-band light-weighted values (right-hand panels) for models of early-type 
galaxi es drawn from 20 realizations of a Coma cluster-size halo. Figure from iTrager & Somervilld 


where A represents the physical parameter (age or metallicity) and mass, and flux, 
are, respectively, the reconstructed mass and flux contributions of the stars in the /-th 
age bin, as returned by the code. When present, young stars are very luminous in 
the optical and, therefore, contribute more to the light-weighted values. This means 
that the light-weighted values of age will be biased towards the youngest stellar 
components. The mass-weighted values will be less biased but they are also more 
uncertain, as the contribution to mass by low-mass faint stars can be very important. 

It is interesting to make a comparison between the SSP-equivalent parameters 
and the averaged ones obtained from th e full s tar formation history. This compar¬ 
ison was made by Trager & Somerviu3 ( 2009h . In their work, the authors derived 
stellar population parameters from synthetic spectra generated by a hierarchical 
galaxy formation model. Figure|^shows the comparison of the SSP-equivalent ages 
and metallicities with the mean values weighted with both light (in the V-band) 
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and mass. As can be seen, the SSP-equivalent ages are always lower than both the 
luminosity- and mass-weighted averages. In particular, SSP-equivalent ages reflect 
more closely the age of the last episode of star formation while luminosity-weighted 
means, although still biased towards the ages of the youngest components, are closer 
to the unweighted mean. On the other hand, SSP-equivalent metallicities and abun¬ 
dance ratios are less severely biased. 

Some studies have used the n on-parametric techniques to analyze the star for¬ 


mation histories of galaxy bulges (I MacArthur et al.[l2()(^ISanchez-Blazquez et al 


2011 : Gonzalez Delgado & et al. . 120141) . The first conclusion from these studies i 


IS 


that the SSP is a ver y bad approxima t ion fo r the star formation history of these 
objects. Furthermore, MacArthur et al. ( 2009h and Sanchez-Blazquez et a DImll) 


found that in a mass-weighted context, all bulges in their sample were predomi¬ 
nantly composed of old stars, independently of their central velocity dispersion. In 
fact, the previously reported trends of age with central velocity dispersion disappear 
when mass-weighted v alues of age are u s ed ins tead of light weighted ones. This can 


be seen in Fig. |9] from MacArthur et alJ ( 2009h . 


The result that bulges are dominated, in mass, by old stars applies to all types of 
bulges (early- and late-type, showing different Sersic index, and with and without 
bars). However, this type of analysis has been performed on a very low number of 
galaxies. It would be desirable to extend this work to a complete sample of bulges 
covering a range of masses and morphological types. Ideally, one would like to 
quantify the mass contribution of the young component and correlate this with other 
properties of the bulges, such as the mass, the environment, the presence of bars, 
and the spiral arm morphology. This would allow one to study the importance of 
secular versus external mechanisms in building up the central bulges. 


galaxies are those performed on data from the CALIFA survey 

/ JPerez <&; et al., 2013: 

Gonzalez Delsado & et al.. 20141 Sanchez-Blazouez & et al. 

20141). None of these 


studie s has yet especially investigated the population in the bulge, but iPerez & et al 


( 2013h analyze the mass assembling history of the central parts of galaxies com¬ 


pared with the rest of the galaxy for galaxies of all morphological types binned in 
mass. They found that galaxies with stellar mass M* > 5 x lO*** Mq have grown 
quickly their inner part, 5-9 Gyr ago, while lower mass galaxies formed their stars 
more slowly. 

One caveat to all these studies performing full spectral fitting is that they use 
stellar population models with chemica l abun dance ratios scaled to solar (e.g. 
Bruzual & Chariot . 20031: Vazdekis et all 2010h . This implies that the models are 


tuned to the specific chemistry and star formation history of our MW. This is because 
the empirical spectral libraries are limited to those stars in the solar neighborhood. 


ratios have been made for SSPs of different ages and metallicities ( 

Trager et al., 

200C 

Thomas et al.. 

2003; Proctor et al.. 2004; Tantalo & Chiosil 2004 

; Lee & Worthev. 

2005 

Annibali et al. 

. 2 OO 7 I Schiavonl 2007) using a semi-empirical approach, the 


calculation of the entire spectral energy distribution is more challenging. How¬ 
ever, in the last few years, full spectrum fitting models have been extended to in- 
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Fig. 9 Average age and metallicity as a function of velocity dispersion for a sample of face-on 
bulges. Black solid squares represent light-weighted values while red symbols are mass-weighted. 
The black and red dott ed lines are linear regressions to the light- and mass-weighted data. From 
iMacArthur et al.l l l2009h . 


elude variation in the elemental abunda nce patterns, using ei t her theoretical stellar 
libraries or semi-empirical ap proaches ( Coelho et al. . 2007t Walcher et all 20091 : 
Conroy & van Dokkum , 2012h . The first studies using full spectral fitting to derive 


chemical abundance ratios of different element s are starting to appear in the liter¬ 
ature, all using samples of early-type galaxies ( Conroy et al. . 20141: Walcher et al 


20091) . This is, however, a challenging task, due to the large number of parameters 


to fit and the possible degeneracies between them. 


5 Spatially resolved stellar populations in bulges 

Most stellar population studies in bulges have been done using the integrated prop¬ 
erties inside a certain aperture. In case of spectroscopic studies, this aperture com¬ 
monly encloses just the very central parts. However, if we want to have a full under- 
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Standing of bulge formation, it is necessary to gain knowledge of the variations of 
the stellar populations with radius. These variations are intimately connected with 
the dynamical processes that led to the formation of these structures, the degree of 
dissipation, and the possible re-arrangement of material. 

Mergers with gas dissipation or monolithic collapse scenarios predic t steep 
metallicity gradients ( Eggen et'd]^ 1962f Larsoni 1^74; Arimoto & Yoshiil 19871) 


The predictions for sec- 


and strong gradients in [a /Fe] JFerreras & Sif 
ularly formed bulges are more complicated. As they formed from redistribution 
of disk stars, the final metallicity gradient will depend on the original gradient 
in the disk and the scale-len gth of the final bulge and also on the disk heating 
(IMoorthv & Holtzmanl 120061) . However, lower metallicity gradients are expected 
compared to those of the first scenario. Observationally, the Galactic bulge, which 
manifests many characteristics of a peanut-shaped bulge, has a clear vertical metal¬ 
licity gradient, such that th e more metal-rich part of the metallicity distribution thins 


out to wards high latitudes ( Minniti et al. , 19951: Zoccali et al. , 2008t Gonzalez et al 


201 ih . This result has long been taken as a signature for a classical bulge in the MW. 
However, recent results have shown that the stars that have been scattered furthest 
from the disk are the oldest stars and, consequently formed from t he least metal- 
enriched fuel ( FreemaiubOOSHMartinez-Valpuesta & Gerhard! l2013h . The buckling 
process may hence establish a nega tive minor-axis meta llicity gradient (which is 
observed in the MW and NGC 4565 ( Proctor et al.l 2000h . 

Several articles have studied the variation of the spectral features with radius 
in bulges ( Moorthv & Holtzmanl 2006: Jablonka et al.l 20071: Morelli et"^ 2008: 


Perez&^anche^ Mazauezl 2011 : lMac Arthur et al.l 20091: Sanchez-Blazquez et al 


201 it Ganda et al. . 2007h and compared them with stellar population models to ob 
tain either SSP-equivalent parameters or mean values, based on the recovery of the 
star formation history. These studies find that most bulges have SSP-equivalent or 
luminosity-weighted negative gradients in the metallicity, almost no gradients in 
age, and slightly positive or null [a/Fe] gradients. Metallicit y gradients in the bulge 
regions are generally steeper than th ose in the disk region (IMoorthv & Holtzman , 
2006t Sanchez-Blazquez et al. . 201 ih . Figure fTOl shows the distributio n of the SSP- 


equiva lent gradients for a sample of bulges taken from the work of iMorelli et al 

(l2008h . 


These values are similar to the ones f ound in elliptical gala xies, although the 
quanti tative comparison is not that clear. .Tablonka et al.l ( 2007h and Morelli et al 


( 20081 ) do not find any diff erence in the magnitud e of the bulge gradients and those 
of elliptical galaxies, while I Williams et al.l(l2012h . on the other hand, found that the 
gradients in boxy bulges are shallower than those in elliptical galaxies at a given a. 

Several studies have also looked for correlations between the gradients and other 
properties of the galaxies, such as the central a, the luminosity or the mass, to check 
if they are related with the potential well of the galaxy, as the central values are. 
Goudfrooii et al. ( 19991) and Proctor et al.l ( 2000t) found that gradients were corre¬ 
lated with luminosity and ce ntral a, respectively, although from very small samples, 
while .lablonka et al.l ( 200^ found no such correlation. How ever, there seems to be 
a trend for which small bulges have lower gradients (see also lMoorthv & Holtzman . 
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Fig. 10 Distribution of the gradients of age (left), metallicity (central), and [a/Fe] enhancement 
(right) for a sample of galaxies. The dashed line represents the median of the distribution and its 
values is also reported. The solid line represents a Gaussian centered on the median value of the 
distribution. Their ff is approximated by the value containing the 68% of objects of the distribution 
and is noted in the inset label . The green and blue a rrow s show the average gra dient found for early- 
type galaxies and bulges by iMehlert et alj l l2003h and Ijablonka et alj <200711 . respectively (figure 
from lMorelli et alJ <2(W8ll l. 


20061: Gonzalez Delgado & et al. . 2014). This, in principle, could be attributed to 
the fact that secularly formed bulges are more common in low mass galaxies. 

The possible differences between the gradients of bulges with and without bars 
have also been explored in a few works, and none of them find any significant one 


( Moorthv & Holtzman. 20o5 Jablon ka et al.L 2007 : Perez & Sanchez-Blazauezl 201 ih 


However, Moorthv & HoltzmanI ( 2006 ) reported that when a positive age gradient 
was present, it was always in barred galaxies, which could indicate that these ob¬ 
jects have more extended star formati on in their centers due to ba r-driven inflow 
of gas. This result agrees with that of iGadotti & dos AniosI (l200lh . who found a 
greater prevalence of null or positive color gradients in barred galaxies than in 
non-barred galaxies, which they interpret as an evidence for gradients being erased 
by bar-driven mixing . Nonetheless, this has not been confirmed in other studies 


(iJablonka et al.L I2007h . The reason for the discrepancies could be, once again, the 


orientation of the samples. The majority of authors agree that positive age gradients 
are normally the consequence of the the presence of central disks or nuclear rings 
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with recent star formation (e.g. lMorelli et alil2008h . As central disks and rings are 
likely more common in barred galaxies, this can explain the differen ces between 
barred and non-barred galaxies found by Moorthv & Holtzmun ( 20061) and also the 
lack of differences reported in the edge-on sample of Jablonka et al. ( 2007h . as these 
flattened central structures will not contribute to the observed light of the bulge in 
these orientations. 

This interpretation of the age gradients as the consequence of the presence of cen¬ 
tral younger structure is supported by the fact that the mass-weighted age gradients 
are, in the majority of cases, much flatter than the luminosity-weighted or the SSP- 
equivalent ones, indicating the the majority of the stars in the bulge are old and share 
a common age, while a small frac tion of stars concentrated in central structures are 
causing the observed radial trends ( MacArthur_etal^ 2009^ Sanchez-Blazquez et al 


201 it Sanchez- Blazque z & et al.L 20141: Gonzalez Delgado & et ah . 2014). 


Interestingly, Jablonka et al.~(l2007h find that the line-strength indices at 1 
were very similar for all the galaxies, independent of their mass or morphological 
type, and the different gradients come from the differences in their central indices. 

An issue in measuring the gradients of age, metallicity and [a/Fe] in bulges 
could be the contamination of their stellar population by the light coming from the 
underlying disk stellar component. This effect is not important in the galaxy center 
but it can have an enormous impact in the stellar population estimates of the external 
parts where the contribution from the disk to the total light is more important. Dif¬ 
ferent authors have tried to qua ntify one way or the other this contamination from 
the disk to the bulg e light (see Jablonka et al. . 1996 : Moorthy & Holtzman . 20061 : 
Morelli et al.[ 20081) and it does not seem to be yery important, but the number of 
tests is small and disk contaminati on is still an is s ue in the measurement of gradi¬ 


ents. Studies of edge-on galaxies (IJablonka et al.L 120071) do not haye this problem 


although they haye the drawback of being blind to flatter components in the center 
of the galaxy. 


6 The bulge-disk connection 


The age distribution of disks is of obyious importance for constraining scenarios 
of disk-bulge formation. Correlations between the disk and bulge colors haye been 
found by seyeral au t hors ( Peletier_&_BalcelB, 1996t ^ Jongl 19961 Bell & de Jongl 


2000 t Carollo et al.L 2001 ; Gadotti & dos Anioa 2001 ) which stands for both early 
and late type spirals. The similarity in color between inner disk and bulge has been 
interpreted as implying similar ages an d metall i cities for these two component s 


mterpreteq as implying similar ages an a metall i cities tor tnese two component s 
and an implicit eyolutionary connection ( de .Tong . 1996L Peletier & BalcelM 1996 ). 
Howeyer, using only colors one cannot directly transform the correlation in colors 
into correlations of age and/or metallicity. 


The radius that contains half of the total luminosity of the bulge. 
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A correlation between the line-str ength index [MgFe]' a t one d isk scale length 
and in the center was also found by iMoorthv & HoltzmanI (l2006h . They interpret 
this as a correlation between the metallicity of bulge and disk as this index is more 
sensitive to variations of metallicity than age and it is almost insensitive to variations 
of[a/Fe]. 

More recently Sanchez-Blazquez & et al.l ( 201-4 performed a full spectral fitting 
analysis of a sam ple of 62 nearly face-on spiral galaxies observed as a part of the IFS 
survey CALIFA ( Sanchez & et al. . 2012h and showed that the slope of the relation 
between the luminosity-weighted age and metallicity and central velocity dispersion 
is similar for the central parts of the bulges and for the disk at ^2.4 scale-lengths, but 
the bulges show higher luminosity-weighted ages and metallicities. FigurefTTIshows 
the comparison of the mean values of age and metallicity (both weighted with mass 
and light) in the center and at 2.5 scale-lengths of the disk. A Spearman rank order 
tests reflect that, while the metallicities of both components are correlated, the same 
does not happen with the ages. 





[Z/H] [at 2.5 rj 



[Z/H] [at 2.5 rj 


Fig. 11 Compari.son of the average ages (left panels) and metallicity (right panels) weighted by 
both mass (top panels) and light (bottom panels) measured in the center and at 2.5 scale-lengths of 
the disk. Different colors indicate if the galaxy is barred (red), non-ba n'ed (blue) or weakly barred 
(orang e). The solid line represents the 1:1 relation. For more details see lSanchez-Blazguez & et"^ 
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Note that the correlation exists for all early and late-type spirals, spherical and 
pseudobulges, contrary to the connectio n between the radial si zes of bulges and 
disks, that only exists for pseudobulges ( Fisher & Drorv . 20081) . Therefore either 
all bulges formed secularly and some had their bar destroyed or other physical 
processes are responsible for this correlation. Several authors have shown that the 
correlation between the disk and bulge sizes can exist in major and minor merger 
remnants. While minor mergers tend to preserve the original bulge-disk coupling of 
the main progenitor, major mergers are capable of rebuilding a bulge-disk coupling 
from the remnants afte r having destroyed the original structures of the progenitors 
( Ouereieta et al. . 20151 and references there in). 

On the other hand, Morelli et al.l ( 20121) rule out significant interplay between 
the bulge and disk components due to the similarity in the stellar population prop¬ 
erties of bulges hosted in galaxies with very different disks (high- and low-surface 
brightness). 


7 Summary 

The study of stellar populations can provide the critical test needed to understand 
the basic mechanism driving bulge formation. However, a clear picture about the 
stellar population in bulges and their possible correlations with other parameters is 
still not obvious. We have collated the main results obtained from the study of colors 
and spectral characteristics and the comparison with stellar population models and 
can thus summarize them: 

• Bulges show a wide range of SSP-equivalent ages and metallicities. There is a 
trend for which more massive bulges have, on average, older stellar populations 
and higher values of metallicity and [a/Fe]. 

• Bulges with high Sersic index ii tend to be old and have high [a/Fe] but it is not 
clear if this trend is due to the existing trend between mass and ii. 

• There are not strong differences in the stellar populations of bulges and elliptical 
galaxies at the same mass, but the details are of this comparison are still not clear 
and conclusions differ between studies. 

• Bulges are, in many cases, composite systems, with disky, boxy/peanut and clas¬ 
sical bulges coexisting in the same galaxy. They usually host different stellar 
populations. Young stars, when present, are located in central disks and rings. 
These young components are not seen in edge-on samples which has lead studies 
differing in the orientation of the sample to obtain different results. 

• In general terms, there are not very clear differences in the stellar population 
properties of bulges with and without bars, neither in the central values nor in the 
variations with radius. If there are differences, those are only present in massive 
galaxies, and they show up as an excess of young populations when compared to 
unbarred galaxies. 

• When the mass-weighted mean values of age are considered, all bulges, indepen¬ 
dent of their mass, seem to be dominated by old stars. The trends between age. 
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metallicity, and mass become much flatter and almost non-existent. This is true 
for bulges with a variety of structural parameters, such as different Sersic indices 
or surface brightness profiles. 

• Bulges show mild, negative, metallicity gradients and almost null [a/Fe] and 
age gradients. The distribution of the slopes is similar to that found in elliptical 
galaxies. At present, there is no agreement about the possible correlation of these 
gradients and other parameters. 

• There is a correlation between the metallicity of the bulge and the disk but not 
between the ages. 


8 Future prospects 


Although considerable progress has been made in the field, there is still much to be 
done, in order to understand the star formation histories and chemical evolution in 
the bulges of spirals and SO galaxies. The first thing we have to deal with is the fact 
that we cannot make a clean distinction between classical, disky, and boxy/peanut 
bulges, as many of them coexist in the same galaxy. We need to quantify the prepon¬ 
derance of each component and correlate this with other properties of the galaxies, 
like the mass, the environment, the presence of bars, type of spiral arms, etc. We 
need also to understand the physical mechanism that formed each components. For 
example, numerical simulations have shown that bulges with structural character¬ 
istics of pseudobulges can be formed, not only secularly, but also quickly, at high 
redshift, via a combination of non-axisymmetric disk instabilities and tidal interac¬ 
tions or mergers. 

Most advances, therefore, will come from the use of new techniques to derive star 
formation histories that allow us to distingu ish different episodes of s tar formation. 
The future is promising; Ocvirk et al. ( 2008 ) and Coccato et al. ( 201 ih show the fea¬ 
sibility of separating different components, not only in t erms of stellar popu lations, 
but also kinematically. This is illustrated in Fig. [12] from lOcvirk et al.l (l2008h . where 
the age-velocity distribution for the bulge is shown in two positions (in the center 
and outside the bulge dominated region of the galaxy). In this work, the authors 
model the observed spectrum as a sum of 40 components with different ages but, 
contrary to the traditional assumption that all components share the same line-of- 
sight velocity distribution (LOSVD), each component was allowed to have its own 
non-parametric (not necessarily Gaussian) LOSVD. These authors reconstructed the 
age-velocity distribution of two bulge regions of the Sbc galaxy NGC 4030 and were 
able to separate two components, one with a relatively young stellar population (^2 
Gyr) and <7 ^ 100 km/s and an even younger ('^ 500 Myr), and kinematically colder 
(O’ ~ 30 km/s). Its location (outside the bulge-dominated region of the galaxy), kine¬ 
matics, and relatively young age suggest that this is a young internal disk. The null 
mean velocity of the structure is expected for the minor axis of a stellar disk. This 
study demonstrates the feasibility of separating and measuring the age and kinemat¬ 
ics of superimposed galactic components from an integrated light spectrum. This 
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type of analysis will be very useful in characterizing the properties of the different 
components in composite bulges and isolating the properties of their stellar popula¬ 
tions. 


age-velocity distribution age-velocity distribution 




Fig. 12 Age-velocity map reconstruction for the center of the bulge spectrum of NGC 4030 (left) 
and for the bulge spectrum 3 arcsec from the ce nter (right panel), wh ere a disk component appears 
on top of the dynamically hot component. From lOcvirk et all 1200 8h . 


Furthermore, we have mentioned throughout the chapter several results that need 
to be confirmed with larger samples. Spectroscopic studies of a large sample of 
bulges with different characteristics are clearly needed to advance in our under¬ 
standing of bulge stellar populations. Ideally this sample will be observed with an 
integral field unit -to be able to separate morphologically the different subcompo¬ 
nents - and it will contain galaxies in different en vironments. Indeed, t he influence 
of environment has not been studied thoroughly. Morelli et al. ( 20081) analyzed a 
sample of bulges in the Co ma cluster but t hey d id not compare their results with 
bulges in other environments. lPeletier et al.l ( 19991) did not find any difference in the 
colors of galaxies in different environments and concluded that this parameter does 
not have a strong influence in shaping the star formation in bulges. Nonetheless, a 
more comprehensive study remains to be done. A detailed study of the properties of 
bulges with environment might help us to distinguish between external and internal 
processes for the formation of bulges. 


We expect that instruments that are already operating, such as MUSE (iBacon & et al 


20K)h . with superb spatial resolution, will help us to resolve central components. 


such as nuclear disks or bars. These data will allow the coupling of stellar popula¬ 
tion and kinematical properties with the morphological characteristics of the galax¬ 
ies, improving our understanding of bulges in a way it has not been possible before. 
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